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SECTION  I 


INTRODUCTION 


During  the  development  of  a  method  for  unsteady  linearized 
subsonic  flows  the  question  arose,  how  the  vorticity  within  the 
wake  is  distributed  during  the  first  stage  of  the  wake  formation. 
This  problem  is  treated  here  for  both  compressible  and 
incompressible  flows. 


The  incompressible  flow  has  been  explored  without  the 
restriction  to  the  initial  stage  by  Herbert  Wagner  (Reference  1). 
The  present  study  uses  Wagner's  concepts.  With  the  restriction  to 
the  initial  stage  one  obtains  rather  specific  results. 


The  nature  of  the  problem  changes  if  one  takes  compressibil¬ 
ity  into  account.  In  incompressible  flow  the  velocity  of  sound  is 
infinite;  the  flow  field  adjusts  itself  immediately  to  changing 
boundary  conditions  and  a  changing  vortex  distribution,  even  if  the 
changes  occur  within  short  times.  In  compressible  flows  perturba¬ 
tions  travel  during  a  finite  time  only  over  a  finite  distance.  In 
the  beginning  stage  the  wake  formation  has  an  effect  on  only  a 
small  part  of  the  flow  field  in  the  vicinity  of  the  trailing  edge, 
but  within  this  part  the  unsteady  effects  cannot  be  neglected. 
Nevertheless,  the  basic  results  are  the  same  in  the  two  cases. 


In  both  cases  one  obtains,  conceptually,  the  actual  flow 
field  by  the  superposition  of  two  types  of  solutions  of  the  partial 
differential  equation.  One  of  these  solutions  satisfies  the  upwash 
conditions  at  the  wing,  but  does  not  allow  for  wake  vortices.  The 
second  solution  is  obtained  by  a  superposition  of  flow  fields  each 
with  zero  upwash  at  the  wing,  they  are  generated  by  one  vortex  shed 
from  the  leading  edge  at  a  certain  time  and  traveling  from  then  on 
downstream  with  the  flow  velocity.  In  an  incompressible  unsteady 
flow  t.he  potential  is,  of  course,  time  dependent,  but  the  time 
dependence  does  not  appear  in  the  potential  equation  (given  by  the 
hap! ace  equation).  Here  these  flow  fields  can  be  represented  by 
explicit  formulae.  In  the  compressible  case,  the  potential 
•  quation  can  be  brought  into  similarity  form,  so  that  instead  of 
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v  *.  *.  -  v  •  •_  *  •  v  \ 


three  independent  variables  (two  space  variables  and  time)  one 
deals  only  with  two  independent  variables.  But  the  resulting 
partial  differential  equation  is  rather  complicated.  Fortunately, 
only  a  very  limited  amount  of  information  is  needed  to  find  the 
vortex  distribution  within  the  wake.  A  general  discussion  is 
sufficient;  only  one  constant  (dependent  upon  the  Mach  number) 
remains  undetermined.  In  the  incompressible  case  these  data  are 
expressed  in  explicit  formulae.  The  availability  of  such  formulae 
allows  one  to  give  specific  information  about  the  resulting  flow 
field.  This  facet  is  derived  in  an  Appendix  B. 


SECTION  II 


INCOMPRESSIBLE  FLOWS,  BASIC  CONCEPTS 

Considered  is  a  two-dimensional  nonsteady  incompressible 
flow.  Let  x  and  y  be  a  system  of  Cartesian  coordinates  in  which 
the  x-axis  has  the  free  stream  direction,  and  let  t  be  the  time.  A 
thin  profile  is  subject  to  small  time  dependent  motions  or 
deformations  in  a  flow  with  the  velocity  U.  The  wake  vortices  move 
with  the  local  velocity.  The  perturbation  of  the  original  parallel 
flow  due  to  the  wing  are  considered  as  small.  The  vortices, 
therefore,  move  with  the  free  stream  velocity  U.  The  airfoil  is 
replaced  by  a  slit  extending  along  part  of  the  x-axis.  An  upwash 
will  be  imposed  along  this  slit  that  depends  upon  x;  with  respect 
to  time  it  is  given  by  a  step  function.  With  results  for  such  an 
upwash,  it  is  then  possible  to  generate  the  response  to  an  upwash 
which  has  the  same  x  dependence  and  arbitrary  time  dependence. 
Without  loss  of  generality,  one  can  assume  that  the  step  occurs  at 
time  zero. 

The  perturbation  potential  depends  on  x,y,  and  t.  But  the 
potential  equation 


0 


does  not  contain  the  time  explicity.  We  do  not  allow  the  airfoil 
t. hickness  to  change  with  time.  The  upwash  is,  therefore,  the  same 
on  the  upper  and  lower  sides  of  the  airfoil.  It  follows,  that  the 
upwash  4>  is  symmetric  with  respect  to  the  x-axis.  The  potential, 
its  derivatives  with  respect  to  x  and  t,  and  the  pressure  are  then 
antisymmetric.  There  will  be,  of  course,  a  jump  of  the  potential 
between  the  upper  and  lower  sides  of  the  wing  and  of  the  wake.  At 
the  wake  the  pressure  is  continuous,  and  therefore,  because  of  its 
antisymmetry,  zero. 


Because  of  the  assumption  of  small  perturbations,  the 
perturbation  pressure  is  given  by 


Since  at  the  wake  the  pressure  is  zero,  it  follows  that  the 
potential  at  the  upper  side  of  the  wake  has  the  form 

4>(  t ,  x,  0+)  =  f(t  -  (x  -  xtr)/U)  (?) 

where  x^r  is  the  x  coordinate  of  the  trailing  edge.  (0+  means 
that  one  approaches  y  =  0  from  above.)  For  the  lower  side  one  has 
the  same  expression  with  the  opposite  sign. 

The  perturbation  field  is  decomposed  into  two  parts.  A 
steady  part  is  given  by  the  circulation-free  flow  determined  by  the 
upwash  condition  (after  the  step  has  occurred)  at  the  wing.  In 
such  a  flow  one  finds  adjacent  to  the  trailing  edge  an  infinite 
pressure  at  points  of  the  wing  and  an  infinite  upwash  at  points  of 
the  wake.  Superimposed  to  this  steady  flow  is  one  caused  by  the 
time  dependent  vortex  distribution  within  the  wake  (as  it  exists  at 
the  current  time).  At  the  wing  zero  upwash  is  prescribed.  Outside 
the  wake  and  the  wing  this  flow  field  satisfies  the  Laplace 
equation.  This  second  field  gives  singularities  at  the  trailing 
edge  of  the  same  kind  as  the  steady  field.  The  potential  at  the 
upper  side  of  the  wake  must  have  the  form  of  Eq .  (2).  The  function 
f  must  be  chosen  so  that  the  trailing  edge  singularities  cancel 
those  of  the  ci rcul ati on- f ree  flow. 

The  potential  of  the  circulation-free  flow  can  be  developed 
with  respect  to  the  distance  from  the  trailing  edge.  Let  the 
origin  of  the  x,y-system  lie  at  the  trailing  edge  and  let 

z  =  x  +  i  y 

Then  the  lowest  order  terms  of  the  development  of  the  perturbation 
potential  are  given  by 


with 


1/2  3/2 
al  /2  Z  +  31Z  +  a3/2  2 


Since  for  this  flow  the  boundary  conditions  are  independent  of 
time,  the  coefficients  "a"  are  independent  of  time.  No  constant 


term  occurs  because  at  the  location  of  the  wake  there  is  no 


potential  jump  in  the  circulation-free  flow.  Then 


=  Im(dft/dz) 


<j>  =  Re(dQ/dz) 


Specifically  within  the  wake,  i.e.,for  z  =  x  >  0 


<j>  =0 

x 


-1/?  1  /p 

>y  =  (1/2)a1/2  x  ♦  a1  +  (3/2)a3/2  x1 


and  at  the  upper  side  of  the  wing  z  =  x  <  0  (z  =  - } x | ) 


- 1  /? 

<Dx  =  -(1/2)a1/2|x|  +  (3/2)a3/2 


<$>  =  a.  +  2a„x  + 

y  1  2 


.''•poci  fi  0  examples  can  be  readily  obtained,  for  instance,  from  the 
formula*-'  derived  in  Appendix  E  of  Reference  2. 


Notice  that  the  term  i  n  <p  and 

x 


with  the  factor  lx 


-1  /2 


-ire;  with  lh°  same  coefficient.  The  terms  in  <|>  at  the  wing  do 


not  contain  fractional  powers.  (They  are  determined  by  the 
boundary  conditions  for  the  upwash.)  The  wake  development  i; 


solely  determined  by  the  coefficient  a 


Formulae  for  a  single  vortex  in  the  wake  in  the  presence  of  a 
fir. do  wing  are  deri.od  in  Appendix  A.  These  are  the  formulae  on 
wnioh  the  work  of  Wagner  is  based.  In  the  present  discussion  where 


rent  riot,  ourselves  from  the  outset  to  small  times  and  non- 
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sequently  to  small  distances  from  the  trailing  edge  one  obtains 
simplified  formulae,  because  then  the  wing  chord  is  very  large  in 
comparison  to  the  distances  under  consideration.  Therefore,  we 
assume  for  the  evaluation  of  the  wake  potential  that  the  wing 
extends  along  the  x-axis  from  -»  to  zero.  The  formulae  so  obtained 
could  alternatively  be  obtained  by  developing  the  complete 
expression  of  the  appendix  under  the  assumption  that  the  distances 
from  the  trailing  edge  are  small.  To  obtain  the  potential  for  a 
singlQ  vortex  in  the  presence  of  such  an  infinite  slit  we  first 
consider  in  an  x1  y1  plane  (Figure  1) 


♦  (  z  1  ) 


I m  log 


a 

a 


a  >  0  real 


with 


Z1  "  X1  *  iy1 

This  is  the  flow  field  with  two  vortices  of  opposite  sign  at  the 
points  x1  =  a  and  x1  =  -a.  The  velocity  component  normal  to  the 
y^axis  is  zero.  The  y^axis  is  mapped  into  a  slit  from  -■»  to  zero 
along  the  x-axis  by  setting 


f,1  /?  f,  >  0  real 


then  one  obtains 

z'/2  -  r1/?  i/2  U2 

't>  '  Log  - =  I m[  1  og ( z  -  Fj  -  2  log(z  _  +  F,  (3) 

z  +  F, 

I'hi  f‘Xpr--ss  i  on  has  one  logarithmic  singularity  at  z  =  F, .  Other 

1  /') 

i  ngul  nri  ties  occur  at  the  trailing  edge  because  of  the  power  z  . 
■  o  !  t  i  i  nr, 
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THE  FORM  OF  THE  WAKE  POTENTIAL  IN  INCOMPRESSIBLE  FLOW 


We  express  the  potential  due  to  the  wake  by  a  linear 
combination  of  expressions  (3).  Let  f(£)  be  some  function  of  the 
umbral  variable  5.  This  potential  is  then  given  by 


rUt  1/2  _  1/2 

*wake(t’x'y)  '  l  fU,Im  log  T72 - 7T72  d( 

0  z  +  £ 


'  I  1 /p  1  /? 

fU)Im  log  (z  -  Od£  -  2  j  f(C)Im  log(z  +  i  )d£ 


In  the  limits  of  integration  we  have  taken  into  account  that  the 
wake  extends  from  0  to  Ut.  The  second  term  in  the  last  equation  is 
analytic.  The  relation  between  the  function  f  and  the  jump  of  the 


potential  at  a  point  xq  between  the  upper  and  lower  side  of  the 


wake  can  be  determined  in  the  following  manner.  The  circulation 
integral  4*  grad  4>-ds  for  a  single  vortex  Im  log(z  -  £)  is  2ir. 


Consider  now  a  path  which  starts  and  ends  at  the  same  point  xq  of 


the  wake  and  which  encloses  the  wake  downstream  from  the  point  x 


(the  wake  ends  at  a  finite  distance  namely  x  =  Ut).  (See  Figure 
2.)  The  second  term  on  the  right  in  Eq.  (4)  is  regular  at  and  in 
the  region  within  this  path  and  therefore  gives  no  contribution. 


One  obtains 


grad  <t>  ( t ,  x  ,  y  )  .  ds  =  <j>(t,xQf0  )  -  <J>  ( t ,  xQ  ,  0  )  =  2  tt  J  x  fU)dE( 

o 


But  <(>  is  antisymmetric.  Hence, 


4>(  t,  xQ,0  ) 


f  U)d£ 


f(xo)  '  -  1  H<t’xo,0*) 
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Figure  2 


Wing  and  wake  at  time  t.  Path  of  integration  for  the 

determination  of  the  potential  difference  between  the 

upper  and  the  lower  sides  of  the  wake  at  a  point  xq. 

The  path  starts  and  ends  at  xQ.  One  proceeds  around  the 

vortices  enclosed  between  £  «  x  and  £  ■  Ut  in  the 

*  o 

counterclockwise  sense. 


mm. 


The  form  of  the  potential  within  the  wake  has  been  found  in  Eq. 
(2).  The  derivative  3<J>/9x  has,  of  course,  the  same  form.  The 
function  f(£)  therefore  has  the  form  g(t  -  (£/ U)),  and  one  obtains 


■  I 


ut 


g(t  -  U/U))Im  log^-72 - ^72  d£ 

z 


1  /2 


I 
+  5 


1  /2 


(6) 


The  singularity  in  <j>  which  arises  in  this  expression  at 
points  of  the  wake  adjacent  to  the  trailing  edge  must  cancel  the 
corresponding  singularity  in  the  circulation-free  flow  around  the 
wing.  This  is  the  condition  which  determines  the  form  of  the 
function  g.  We  form 


9l^wake 

ay 


ut 

g(t 

0 


U/U))Re 


-1  /2 


z 


For  z  =  x  >  0  one  obtains 


d^wake 
ay  " 


U/ U)) 


In  order  for  this  expression  to  cancel  the  corresponding 
singularity  in  the  circulation-free  flow  one  must  have 

fut  r1/2 

1  =  x-0  '  g(t  '  U/u)  d£  =  const 


where  the  right-hand  side  does  not  depend  upon  the  time.  Let 


S  =  q  Ut 

x  =  x  Ut 


Then  one  obtains 


1 


1  im  ,,1/2  |  .1/2 
x  +  0  0 


1  /2 


U”fcJ  X?/Z  g(t(1  -  q))^ -  dq 


x  -  q 


1  1 


and  one  postulates 


.  .  f  ,  /0  1/2 

J  t  2g(t(1  -  q) )  dq  -  const 

x+0  0  x  -  q 

1  /2 

This  gives  the  requirement  that  t  g ( t ( 1  -  q))  be  solely  a 

function  of  q. 

g(t(1  -  q))  =  f(q)t"1/2 

The  argument  of  g  is  t(1  -  q).  It  follows  that  f(q)  has  the  form 
-1/2 

const(1  -  q)  .  Therefore, 

g(t(1  -  q))  -  t  1/2(1  -  q)  1/2  const 

or,  after  one  returns  to  the  original  coordinate  with  a  different 
choice  of  the  constant, 

-I/O 

g  -  const(Ut  -  £)  (7) 

Thus , 

f Ut  /?  1/2  _  1/2 

<J>wake(t,x,y)  =  const  J  (ut  “  €)  1  Im  1°8Z1  /2  y/2  (8) 

0  z  +  £ 

Eq .  (8)  is  the  crucial  result  for  the  incompressible  flow. 


SECTION  IV 

THE  WAKE  IN  LINEARIZED  COMPRESSIBLE  FLOW 


As  in  the  incompressible  flow,  one  suddenly  imposes  at  time 
zero  the  boundary  condition  of  constant  upwash  at  the  wing.  If  the 
wing  extended  along  the  entire  x-axis,  one  would  obtain  a  compres¬ 
sion  wave  and  an  expansion  wave  respectively  on  the  upper  and  lower 
sides  which  at  a  (positive)  time  t  generate  fields  of  positive  and 
negative  perturbation  pressure  and  constant  upwash  out  to  a 
distance  a  t.  (See  Figure  3.) 

Perturbations  float  downstream  with  the  free  stream  velocity 
U  and  expand  with  the  velocity  of  sound  a.  If  the  plate  ends  at 
x  *  0 ,  the  above  flow  field  will  terminate  at  a  circle  around  the 
point  x  =  Ut  with  radius  a  t.  We  consider  times  which  are 
sufficiently  small,  so  that  perturbations  coming  from  the  leading 
edge  do  not  affect  the  region  in  the  vicinity  of  the  trailing 
edge.  (See  Figure  4.) 

Assume  first  that  no  vortices  are  shed  from  the  trailing 
edge.  Within  the  circle  described  above  one  then  obtains  a  flow 
field  caused  by  the  pressure  difference  between  the  upper  and  lower 
side.  This  flow  field  will  have  similarity  form  with  respect  to 
time,  that  is  the  velocities  will  depend  only  upon  x/t  and  y/ t .  At 
the  portion  of  the  plate  within  the  circle  one  will  have  the 
desired  upwash.  At  the  portion  1,  2,  3  of  the  circle  (which  moves 
with  time)  the  solution  within  the  circle  must  match  the  solution 
for  the  infinite  plate,  at  the  remaining  portion  of  the  circle  it 
must  match  the  undisturbed  flow.  One  expects  that  at  points  of  the 
x-axis  downstream  of  the  trailing  edge  and  adjacent  to  it  this  flow 
has  a  singularity  in  the  upwash  of  the  same  character  as  in  an 
incompressible  flow.  This  will  be  discussed  in  some  detail. 

Superimposed  to  this  field  is  another  one  due  to  the  vortices 
shed  from  the  trailing  edge.  A  single  vortex  which  leaves  the 
trailing  edge  at  time  t  moves  downstream  with  the  velocity  U.  It 
generates  a  field  which  also  has  similarity  character;  it  depends 
solely  upon  x/(t  -  t)  and  y/(t  -  t).  The  intensity  of  the  vortices 


1 - ? - ? - TT 

CONSTANT  UPWASH  at 

I  t  I _ il 

1 - 1 - 1 - TT 

CONSTANT  UPWASH  ot 

1 _ 1 _ 1 _ il 


COMPRESSION  WAVE 


INFINITE  PLATE 


EXPANSION  WAVE 


Figure  3.  Perturbation  generated  in  a  parallel  flow  with 
velocity  U  by  an  infinite  plate  (double  line) 
suddenly  set  in  motion  in  the  y-direction  with 
a  constant  velocity.  One  obtains  a  constant 
upwash  for  |y|  <  at  and  no  upwash  at  |y|  >  at. 
At  y  =  at  and  y  =  -at  one  has,  respectively, 
a  compression  and  a  rarefaction  wave. 
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shod  in  this  manner  will  depend  upon  t.  The  function  which  gives 
the  intensity  must  be  determined  in  such  a  manner  that  the 
singularity  generated  in  the  circulation-free  flow  is  cancelled  at 
all  times  by  the  flow  field  generated  by  the  vortices  shed  form  the 
trailing  edge.  In  the  following  we  shall  develop  these  ideas  in 
detail . 

After  a  linearization  for  the  vicinity  of  a  parallel  flow  the 
perturbation  potential  in  the  two-dimensional  case  satisfies  the 
equat i on . 

(a2  -  U2)»xx  ♦  a2*yy  -  2U„xt  -  «tt  -  0  (9) 

The  potential  of  a  vortex  moving  downstream  with  the  velocity  U  (in 
other  words,  with  zero  velocity  with  respect  to  the  surrounding 
fluid)  is  gi ven  by 

4)  *  Im£l 

n  -  1 og ( x  +  iy  -  Ut)  (10) 

The  real  part  in  Eq .  (10)  would  give  a  moving  source.  This 
expression  should  satisfy  Eq .  (9).  Let 

z  -  x  ♦  i  y 


Similarity  solutions  are  obtained  by  introducing  independent 


variables 


5  =  x/t  ,  n  =  y/t  ,  t  =  t 


(11  ) 


4>(  x ,  y ,  t )  =  t  i|>(  £.  n,  t  ) 


(12) 


If  41  actually  depends  upon  the  third  variable  t,  this  is  merely  a 
coordinate  transformation.  A  similarity  solution  arises  if  the 
function  4<  is  independent  of  t.  To  describe  the  circulation-free 
flow,  we  must  set  a  =  1.  It  gives  for  fixed  5  and  n  velocities 
which  are  independent  of  time  (t  or  t).  This  is  in  accordance  with 
the  boundary  condition  of  constant  upwash.  In  the  velocity  field 
due  to  a  single  vortex  moving  with  the  velocity  U  one  must  set  a  = 

0  for  then  the  circulation  around  this  vortex  (the  jump  of  <p)  is 
constant  if  one  travels  around  the  vortex  along  some  path  returning 
to  the  same  point. 


One  obtains  from  Eq .  (12) 


=  T  t|>. 


=  T  4\ 


TU  { (  a-1  )  ipr  -  54v  +  Ti  ) 


t  {a(a-1)4>  +  2(a-1)(-£\p  -  rnj;^)  +  2aTt|>T 


+  ^*55  +  2f’n^5n  +  n\n  *  '  2nT,|;nT  +  T\T: 


Then  one  obtains  from  Eq.  (9)  for  a  =  0 


-  b'2  ♦  25U  -  £2)  +  ^--(a2  -  n2)  +  \pF  n(2Un  -  25n) 


( 1  3) 


vP(2U  -  25)  ♦  ip  ( -  2  n )  ♦  xip  (  — 2U  +  25)  +  rip  (2n)  -  t  \p  =  0 
n  t,  t  n  t  r  t 


Futhermore,  f or  a  =  1  and  independent  of  x 


4^(a2  -  U2  ♦  2£U  -  £2)  +  ^nn(  a2  -  n2) 


+  ^Cn(2Un  -  2£n)  =  0 


(14) 


The  last  two  equations  have  the  same  principal  parts  (coefficients 
of  the  highest  derivatives). 

The  behavior  of  particle  solutions  at  the  trailing  edge  can 
be  studied  by  a  development  with  respect  to  £  and  n.  As  mentioned 
above,  <p  is  antisymmetric.  The  form  of  the  lowest  terms  is 
obtained  by  setting  E,  =  0  and  n  =  0  in  the  coefficients  of  Eqs. 

(13)  and  (14).  Suppressing  the  time  dependence  in  4)  one  obtains 

4>  _(a2  -  U2)  ♦  <p  a2  -  0  for  o  -  1  (15) 

£.£,  nn 

^(a2  -  U2)  ♦  t|>nna2  +  2U^  =  0  for  a  =  0  ( 1 6 ) 

The  term  Uij>  in  the  second  equation  already  gives  a  contribution  of 
higher  order  in  ijj.  In  essence  one  deals  with  the  Laplace  equation. 
Setting 

q  =  U  ♦  i(1  -  M2 ) ^  /2n)  ,  M  =  U/a 


'j  n  o  obtains  for  the  lowest  order  terms  in  the  development  in  Eq . 

(15) 

,  T  1  /2 

=  a1  /2  I  me  +  a  Im; 

The  upwash  in  the  case  a  =  1  is  then  given  by 


d  4> 

d  y  S  n 


M2.1/2r,  /0sD  ,  -1/2,  , 

M  )  [(a1/2/2)Re(;  )  ♦  a1  ] 


-1/2 

1 1'"’  wir.g  (r  <  0),  Re(^  )  =  0.  The  coefficient  a^  is  given  by 
‘u"  boundary  condition  for  <j>  prescribed  at  the  plate.  So  far  the 


'•*  ^  ^ 


-»*L-  Wl- 


coefficient  ^  is  unknown.  Its  determination  requires  the 
solution  (in  terms  of  £  and  n)  of  the  boundary  value  problem  for 
the  circulation-free  flow. 

The  flow  due  to  a  single  vortex  moving  in  the  xy-system  with 
a  velocity  U  is  required  to  satisfy  the  boundary  condition  =  0 
at  the  wing.  Its  behavior  in  the  vicinity  of  the  trailing  edge  is 
expected  to  be  given  by 

<t>  =  bW2  lm(  ;  1  ) 

(The  coefficient  corresponding  to  a^  is  zero.) 

To  obtain  an  expression  which  gives  the  moving  vortex 
singularity  we  rewrite  F.q .  (10) 

=  Im[log(*  ♦  y*-  -  U)  +  *°g  i  J  =  Im[log(c,  ♦  in  ~  U)  *  log  ;  j 

( 1  i  ) 

Except  for  the  additive  term  log  t  this  expression  nas  similarity 
form.  The  entire  expression  satisr  :  es  tv; .  .1  3 )  ,  because  the 

original  expression  satisfies  E  3 .  (  >  .  I  n*>  r  *-.i  1  arid  imaginary 
parts  of  the  expression  (17)  subs  1 1  ♦  at  <•  1  into  Eq .  (1  j)  will  give, 
respect  i  vely ,  a  real  and  an  imaginary  expression.  But  Im  1  of:  t  •*  0 . 


T  h erefore 

the 

imaginary  part  s  a  t  i 

sfies  Eq .  ' 

1  iJ  , 
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Details  are  omitted.  Accordingly,  the  expression  describing  a 
moving  vortex  is  given  by 

4)  =  Im  log(E  +  in  -  U)  (18) 

This  expression  fails  to  satisfy  the  boundary  conditions,  zero 
upwash  at  the  wing,  and  zero  perturbation  at  the  circle  around  the 
point  E,  =  0,  n  =  0  with  radius  "a."  Therefore,  a  particular 
solution  is  superimposed  which  corrects  for  this  failure.  It 
depends  on  the  parameter  U/a  =  M.  At  the  trailing  edge  it  will 
have  a  singularity 

b1 /2(M)Im(E  ♦  (1  -  M2)1 72  in) 1 72 

The  value  of  b1/2(M)  is  the  only  information  needed  to  determine 
the  potential  in  the  wake.  Without  detailed  computation  the  value 
of  this  constant  is  not  available. 

We  return  to  the  coordinates  (x,y,t).  The  expression  with 
zero  upwash  at  the  wing  for  a  vortex  generated  at  the  trailing  edge 
at  time  t,  has  a  logarithmic  singularity  at  £  =  U,  n  =  0  given  by 

<j>  =  Im  log( — ~ 7  +  — -  U) 
t  -  t  t  -  t 

A-  the  origin,  this  particular  solution  gives  a  term 


Let  the  intensity  of  the  vortex  shed  at  time  t  be  g(t).  The  time 
dependent  singularity  at  the  trailing  edge  is  then  given  by 


,  r  ,  .  u2,1  / 2  vl  /2 

b(M)’/2  I  g(t)Im  (  1-^-2 - I)  , 

0  N  t  -  t 


k/m»1/2t  ,  M2,1/2  ,1/2  [ 

=  b(M)  Im(x  +  i(1  -  M  )  y)  J 


0  (t  -  t)1  /2 


Now  the  postulate  is  imposed  that  this  expression  cancel  the 
singularity  in  the  circulation-free  flow  which  was  given  by 

a1/2(M)Im(x  +  i(1  -  M2)1/2y) 

This  expression  does  not  depend  upon  time.  The  integral 


I  -*LPt7S 

o  (t  -  t y/d 


must  therefore  be  a  constant,  independent  upon  time.  Let 


1  -  f  "  v 


The  integral  then  transforms  into 


j  v  1 /2[t1 /2g(t(1  -  v) ]dv 


In  order  for  t  g ( t ( 1  -  v))  to  be  independent  of  t,  one  must  have 

.  .  - 1  /2  ,  .  ,-1/2 

g  =  const  t  ( 1  -  v) 


or  after  substitution  of  v 


„  t  -1/2 

g  =  const  t 


This  gives  the  time  history  of  the  vortex  shedding.  Because  of  the 
antisymmetry  of  the  potential  with  respect  to  the  x-axis,  the 
circulation-free  flow  and  the  potential  added  to  the  expression 
( 1 8 )  (to  satisfy  the  boundary  conditions)  give  zero  potential  for 
x  >  0  (for  x  <  0  at  the  wing  there  will,  of  course,  be  a  potential 
jump).  The  potential  as  one  approaches  the  wake  from  above  or 
below  is,  therefore,  solely  given  by  the  superposition  of 
expressions,  Eq.  (18). 

t 

4>(  t ,  x  ,  +0  )  =  J  g(  t )  Im  log  (X  *  iy  -  U\dt 
0  U  -  t  / 

A  clearer  picture  for  a  fixed  time  t  is  obtained  by  introducing 

U  ( t  -  t )  *  x 


Then  one  obtains,  with  a  different  constant, 


Ut 


$( t ,  x , +0)  =  const 


1 


0  (Ut  -  x) 


- . f/2  Im(log(x  +  iy)  -  x)dx 


which  shows  the  vortices  which  appear  at  the  station  x.  The  vortex 


distribution  extends  from  zero  to  Ut,  the  intensity  is  given  by 
—  —1/2 

eonst(Ut  -  x)  .  According  to  the  discussion  given  in  conjunc¬ 


tion  with  the  incompressible  case  ( Eq .  (5)),  one  then  finds 


<|»(  t ,  x  , +0  )  =  2n  const(Ut 


x) 


1  /2 


-1  /2 

The  coefficient  of  the  x  singularity  for  at  x  <  0  and  for  <p 

x  y 


at  x  >  0  in  the  expression  due  to  the  shed  vortices  is  independent 
of  time.  The  constant  in  all  the  expressions  is  chosen  so  that 
they  cancel  the  corresponding  term  in  the  circulation-free  flow. 

Notice  the  (Ut 

downstream  end  of  the  wake.  It  is  caused  by  the  fact  that  the  wake 
ends  rather  abruptly.  The  potential  goes  down  to  zero  as 


-1/2 

x)  singularity  that  arise  at  the 


1  /2 

(Ut  -  x)  .A  sudden  jump  of  the  potential  would  correspond  to  a 
single  vortex,  and  then  <j>y  would  behave  as  (Ut  -  x)  ^  .  Under 
present  conditions  the  transition  is  somewhat  smoother,  and  the 
power  for  is  -1/2  rather  than  -1. 
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APPENDIX  A 


THE  FLOW  FIELD  OF  A  POTENTIAL  VORTEX  IN  THE  VICINITY  OF  A  PLATE 
IN  TWO-DIMENSIONAL  INCOMPRESSIBLE  FLOW 

The  flow  field  is  obtain  by  a  sequence  of  conformal 
mappings.  Figure  A.1. 

z1  *  V  ♦  iy. 

We  begin  (as  in  the  main  text)  with 

z.  -  a 

$ ( z  1  )  -  Im  log  z---a  ■  a  >  0,  real 


The  transformation 

1  ♦  z1  z2  ‘  1 

Z2  "  1  -  Zl  *  Z1  “  z 2  “f 

maps  the  right  half  of  the  z^-plane  into  the  outside  of  the  unit 
circle  in  the  z2~plane.  (The  left  half  is  mapped  into  its 
i nsi de  . ) 


One  obtains 

z2  -  b 

4>,  (z9)  =■  Im  log  - 

z  *  z2  -  b 

wi  th 


The  transformation 

z  -  (1/2)(z2  +  z2_1  )  ,  z2  -  z  ♦  /z2  -  1 

maps  the  outside  of  the  unit  circle  in  the  z2~plane  into  the 
whole  z-plane  with  a  slit  extending  from  -1  +1.  One  obtains 
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For  the  constant  b  in  the  above  formula  we  now  introduce 


b  -  & 


then 


b 


-1 


Eq.  ( A . 1 )  must  be  applied  if  the  extension  of  the  wake  is  not 
in  comparison  to  the  span  of  the  wing.  It  is  the  expression 
underlying  the  work  of  Wagner. 


small 


APPENDIX  B 


DETAILED  EVALUATIONS  FOR  THE  INCOMPRESSIBLE  CASE 

The  constant  occurring  in  Eq.  (8)  (for  incompressible  flow) 
is  found  from  the  requirement  that  the  expression  <j>  .  cancel  the 

corresponding  singularity  in  the  circulation-free  flow  around  the 
profile.  This  requires  the  evaluation  of  the  integral  occurring  in 
this  equation.  Beyond  this  we  are  interested  in  the  properties  of 
the  flow  field  due  to  a  .  ,  in  particular  along  the  x-axis. 

W  a  KG 

There  may  be  a  question  whether  Eq.  (6)  (which  is  more 
general  than  Eq.  (8))  gives  zero  pressure  within  the  wake. 
Therefore,  we  evaluate  4>x  and  <t> t  directly  from  this  equation.  The 
result  to  be  expected  is  found  from  Eq.  (5) 

4>x  -  -irg(t  -  (x/U))  0  <  x  <  Ut  ( B.  1  ) 

Then,  because  the  wake  perturbation  pressure  is  zero,  from  Eq.  (1) 

<|>t  -  ttU  ( g  ( t  -  (x/U))  0  <  x  <  Ut  ( B.  2 ) 

We  restrict  ourselves  to  an  evaluation  of  4>t.  To  avoid  in  Eq.  (6) 
a  differentiation  of  g  with  respect  to  t,  we  introduce,  in  essence, 
instead  of  £  the  argument  of  g  as  a  new  variable. 


Ut  -  £  -  v 
-d£  =  dv 

The  the  limit  £  =*  0  gives  v  -  Ut, 

the  limit  ^  -  Ut  gives  v  -  0. 

One  obtains 


(B.  3) 


wake 


Ut 

I 


g(  v/U  )  Im  log  —j2 


z1/2  -  (Ut  -  v)1/2 


+  (Ut  -  v) 


T77 


d  v 


29 


I »**•  •%  •*,  /-  /.  »*,  *  .  ' 


The  differentiation  of  the  upper  limit  with  respect  to  t  gives 

1  /2 

Ug(  t )  Im  log^y-72 


The  imaginary  part  of  log  1  is  not  one-valued.  Consider  a  point 
z  =  x  >  Ut.  This  is  a  point  outside  of  the  incipient  wake  and 
because  the  antisymmetry  of  the  potential  with  respect  to  the  x- 
axis  <p  =  0.  Therefore,  log  1  must  be  taken  equal  to  zero  at  these 
points  and,  by  continuation,  everywhere  else.  We  are,  therefore, 
left  with  differentiations  under  the  integral  sign.  One  obtains 

“^wake 

~  3 1 

3<^wake 
3  x 


Ut 


-  I 


g  (  v  /  U  )  - 
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For  z  =  x  >  0,  but  different  from  Ut  -  v,  the  imaginary  parts  are 
obviously  zero.  A  contribution  to  <j>t,  therefore,  comes  only  from 
the  immediate  vicinity  of  the  point  v  =  Ut  -  x.  There  one  obtains 
as  dominant  terms 
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a  >  0  small,  b  >  0  small,  a  and  b  real.  Assume  that  z  approaches 
the  point  x  from  within  the  upper  half  plane. 


Then  one  has 
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Hence  in  accordance  with  Eq.  (5) 
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The  procedure  for  -  -- —  is  in  essence  the  same. 


3  x 


For  the  remaining  evaluations  the  specific  form  of  g  ( Eq 
(8))  is  substituted.  Introducing  v  as  above  (Eq.  ( B . 3 ) )  one 
obtai  ns 
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We  make  the  lengths  dimensionless  with  Ut 
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We  set  w  =  q  .  Then 
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can  be  evaluated  in  terms  of  elementary  functions.  The  basic 
formula,  which,  of  course,  can  be  verified  by  differentiation,  is 
round  in  Reference  2  (Eqs.  236,  3c,  and  3d),  x  <  a  is  real. 
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For  real  a  the  first  version  is  practical  for  |ot|  <  a;  it  can  be 
used  also  for  complex  a.  The  expressions  (B.7)  and  (B.8)  can  also 
be  evaluated  by  the  calculus  of  residues.  The  necessary  trans¬ 
formations  are  shown  in  Appendix  C,  but  in  principle  the  calculus 
of  residue  requires  the  same  transformations  as  the  systematic 
derivation  of  Eqs.  (B.9).  In  the  present  case  a  =  1  and  x  =  q. 
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For  z  =  x  >  1  we  set 
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Then  f rom  Eq .  ( B . 7 ) 
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Here  Eq.  (B.lOa)  is  applied.  We  set 

a  =  ±  i  b 

i  2  .  .  2 

1  -  a  =  1  +  b  »  x 
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After  substitution  of  the  limits,  the  terms  left  within  the  braces 
are  log(-ib)  -  log(ib)  =  -in.  Substituting  b  ( Eq .  (B.11))  one 
obtai ns 
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For  0  <  z  =  x  <  1  one  must  delay  setting  z  =  x  because  of  the 
singularity  which  arises  in  the  integrand  at  q  =  1  -  x.  We  write 
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The  singularity  of  the  integrand  expresses  itself  by  the  term 
1 

log(q  -  a)  .  We  let  z  approach  the  point  x  in  the  upper  half  plan. 
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is  then  a  point  of  the  lower  half  plane,  the  variable  of  integration 
q  is  real  and  ranges  from  zero  to  one.  See  Figure  B.l.  A  cut  from 
q  =  a  to  q  =  -i^  needed  to  make  log(q  -  a)  one  valued,  does  not 

intersect  the  path  of  integration.  In  the  limit  e+0  the  modulus  of 
q  -  a  is  0  for  q  =  1  and  +-n  for  q  =  0.  Figure  B.l  therefore 
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Figure  B.l.  Determination  of  the  argument  of  (q  -  a)  and  (0  -a) 
in  the  complex  q  plane.  A  cut  from  the  point  a  to 
negative  imaginary  infinity  does  not  intersect  the 
path  of  the  integration  from  q  =  0  to  q  =  1 .  As  one 
moves  in  an  x,  y  plane,  (y>  0)  to  y  =  0,  a  moves  in 
the  lower  half  plane  toward  the  real  q  axis. 
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In  the  remaining  terms  one  can  replace  immediately  a  by  a(x) 
A  -  x .  One  obtains 
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For  z  =  x  <  0  one  can  immediately  replace  z  by  -|x|.  Then 
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Here  |a|  >  1  and  the  second  version  of  Eq.  (B.10)  is  applied 
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Substituting  the  limits  one  obtains  the  same  expression  as  for 

0  <  x  <  1 
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We  repeat  Eq.  (B.12) 
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APPENDIX  C 


TRANSFORMATIONS  OF  THE  EXPRESSION  I1 ,  EQ.  (B.7) 
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generates  a  rational  integrand  and  maps  the  points  q  -  -1  and 
q  =  +1  into  p  -  0  and  p  *  •,  respectively.  One  obtains 
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In  this  form  the  integral  is  suitable  for  the  evaluation  by  the 
calculus  of  residues.  The  denominator  can  be  rewritten 
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